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ABSTRACT 

I present near-infrared J (1.25^m), H (1.65^m), and K (2.2/im) imaging observa- 
tions of 185 square arcminutes in 21 high galactic latitude fields. These observations 
reach limiting magnitudes of J ~ 21 mag, H ^ mag and K ~ 18.5 mag. The de- 
tection efficiency, photometric accuracy and selection biases as a function of integrated 
object brightness, size, and profile shape are quantified in detail. I evaluate several 
popular methods for measuring the integrated light of faint galaxies and show that only 
aperture magnitudes provide an unbiased measure of the integrated light that is inde- 
pendent of apparent magnitude. These J, and K counts and near-infrared colors are 
in best agreement with passive galaxy formation models with at most a small amount 
of merging (for VLm = 0.3, JIa = 0.7). 

Subject headings: galaxies: evolution - galaxies: photometry - galaxies: formation - 
cosmology: observations 



Introduction 



The surface density of galaxies as a function of integrated brightness, also known as the 
number-magnitude relation, is one of the classic tests of observational cosmology. This test has 
a long, distinguished history as a tool for studying the nature of the universe and the evolution 
of galaxies (e.g. Hubble 1934; Sandage 1961; Tinsley 1977). Surveys to progressively deeper ap- 
parent magnitude limits have measured the number-magnitude relation at a range of wavelengths, 
and this work has been reviewed by Koo &: Kron (1992) and Ellis (1997). The observed number- 
magnitude relations are compared to the predictions of galaxy evolution models, such as those by 



^Based on observations obtained at MDM Observatory, operated by Columbia University, Dartmouth College, the 
University of Michigan, and the Ohio State University 

^Current Address; Carnegie Observatories, 813 Santa Barbara St., Pasadena, CA 91101, martini@ociw.edu 



- 2 - 



Yoshii &; Takahara (1988), Guiderdoni 8z Rocca-Volmerange (1990), and Gardner (1998) in order 
to determine cosmological parameters and study galaxy formation and evolution. 

Near-infrared (1 - 2.5jLtm, hereafter NIR) galaxy counts are a valuable addition to studies of 
galaxy formation and evolution with the number-magnitude relation. The spectral energy distri- 
butions (SEDs) of most galaxies are relatively constant as a function of wavelength in the NIR 
because this wavelength regime is dominated by light from old, evolved stars. The shape and uni- 
formity of galaxy SEDs in the NIR results in relatively small ^-corrections for galaxies up to 2; ~ 1 
and therefore the NIR luminosity is representative of the total stellar luminosity. A recent burst of 
star formation, in contrast, can significantly increase the brightness of a galaxy at visible and UV 
wavelengths. This may cause a galaxy with a large amount of current star formation to appear as 
bright as a more quiescent galaxy with a much larger stellar mass. The presence of a great deal 
of dust will also attenuate the brightness of galaxies at these wavelengths and further complicate 
the relation between UV or visible luminosity and total stellar luminosity. NIR number counts to 
a limit of ~ 20 mag, which are dominated by galaxies at z < 1, provide a less uncertain means of 
evaluating models of galaxy formation. The old, evolved stars that dominate the NIR light are also 
good tracers of the total stellar mass and therefore NIR measurements are an excellent method for 
direct study of the stellar mass evolution in galaxies in the context of hierarchical galaxy formation 
(Brinchmann & Ellis 2000). 

NIR number counts to date have mostly been obtained at K (Gardner et al. 1993; Cowie et 
al. 1994; Djorgovski et al. 1995; McLeod et al. 1995; Gardner et al. 1996; Moustakas et al. 1997) as 
this is the longest-wavelength atmospheric transmission window with sufficient sensitivity to detect 
high-redshift galaxies in a reasonable time. Faint counts have also been obtained at H (Teplitz et 
al. 1998; Yan et al. 1998) with NICMOS on the Hubble Space Telescope (HST), though given its 
small field of view, there are few bright galaxies in this sample. Some observations have also been 
obtained at both J and K (Bcrshady et al. 1998; Teplitz et al. 1999; Vaisanen et al. 2000) and 
these studies have used the mean NIR colors as a function of apparent magnitude as an additional 
constraint for galaxy evolution models. 

I have obtained J, H, and K observations of 185 square arcminutes of high galactic latitude 
fields which extend to J ~ 21, ~ 20, and K ~ 18.5. The median redshift to K = 18.5 is 2; ~ 0.5 
in the K-selected Hawaii Redshift Survey (Cowie et al. 1996). These observations are made up of 
21 subfields of the Deep Multicolor Survey (DMS) fields observed by HaU et al. (1996a). The DMS 
covers a total area of 0.83 square degrees in six filters: UBVRI-j^I^q. The DMS is comprised of 6 
different fields and to date this rich dataset has been used to study the quasar luminosity function 
(Hall et al. 1996b; Kennefick et al. 1997), the evolution of blue galaxies (Liu et al. 1998), and the M 
dwarf distribution, mass, and luminosity function (Martini k, Osmer 1998). The selection criteria 
for the DMS fields were relatively low galactic foreground extinction and avoidance of very bright 
stars. Five of the fields are equatorial to facilitate follow-up spectroscopy from both hemispheres. 

In §2 I describe the NIR observations, while in §3 I discuss the image processing, photometric 
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solution, and catalogs. I discuss the detection efficiencies in §4, galaxy photometry in §5, and 
star /galaxy separation in §6. The galaxy number counts and NIR colors of these galaxies are 
presented in §§7 and 8. In a separate paper (Martini 2000) I combine these NIR observations with 
the UEVRIf^Ise CCD data to compute photometric redshifts for this NIR sample to study the 
visible-NIR colors of these galaxies, particularly the Extremely Red Objects. 

2. Observations 

A total of 21 randomly-selected subfields of the DMS were observed through the NIR J, 
H, and K filters. These images were obtained with TIFKAM^ (Pogge et al. 1998), a NIR im- 
ager/spectrograph with a 512x1024 InSb detector on the 2.4m Hiltner telescope of the MDM 
Observatory. The f/7 camera was used for all of the observations; this camera gives a plate scale 
of 0.3" pix^^ on the Hiltner telescope. Table 1 is a log of the observations, listing the date each 
field was observed in each filter and the seeing measured in that field/filter combination. Each 
of the 21 fields was observed for a total of 60 minutes on-source integration time per filter. The 
observations were obtained in a six-point dither pattern with a 10" - 15" offset between positions. 
Individual exposures in the dither pattern were one to two minutes and comprised a combination of 
integration time and coadds such that the sky level in the frames was less than 10000 counts. Tests 
of the detector linearity showed that the array is less than 1% nonlinear at this count level. Dark 
frames were taken nearly every night with the same combination of exposure time and number of 
coadds as the science frames. 

Standard stars from Persson et al. (1998) were observed each night at a range of airmass. A 
series of red stars was also observed to determine the transformation of the instrumental magnitude 
system to the CTIO/CIT photometric system (Elias et al. 1982). Red stars were observed at least 
one night per observing run as there were several changes to the filters over the course of the 
observations for this program. These changes are discussed below in §3.2. The standard stars were 
observed in a five or six point dither pattern with 15" offsets, though with only 15 - 30 seconds 
on-source integration per dither position. 



^TIFKAM: The Instrument Formerly Known As MOSAIC. This instrument is also known as ONIS at Kitt Peak. 



-4- 



3. Image Processing 
3.1. Survey Fields 

All of the data were reduced with IRAF^ and many of the data processing steps took advantage 
of the PHIIRS^ package, a collection of IRAF routines compiled by Pat Hall. The first step in the 
data processing was subtraction of a dark frame from each image. The individual images were 
then flat-fielded and sky-subtracted using "running" fiat and sky frames. The running fiats were 
created by averaging the six neighboring images using a percentile clipping algorithm to remove 
bright objects. This image was then normalized to unity to create a skyflat for each image. The 
running sky frames were similarly created by averaging the six neighboring frames. At this stage 
the images were inspected and images with anomalously high noise or serious tracking/guiding 
errors (generally two to three per field) were removed from further processing. The good images 
were then shifted and coadded together using an offset table of integer shifts to create a final mosaic 
image. This image was used to identify objects and, together with the offset table, was used to 
make an object mask for each of the individual frames. The flat field and sky subtraction steps 
were then repeated using these individual object masks. 

TIFKAM has a noise pattern in the lower half of the array due to vibration of the mechanical 
cryocooler. This pattern is a few counts in amplitude and runs along rows in the form of a sine 
wave with a period six times the width of the array. Because the period is much longer than the 
array width and the amplitude is less than the sky noise, simply subtracting the average of each 
row, after masking out objects, removes most of this feature. After the cryocooler pattern was 
subtracted, the individual images were shifted and added into the final mosaic. This region was 
then trimmed to only include the area of the sky present in more than 83% of the individual frames. 
As the offsets were generally small and consistent, this step produced a final mosaic of relatively 
constant noise with minimal loss of field. Several of the fields in Table 1 were initially observed 
under nonphotometric conditions. Calibration images of these fields were obtained on photometric 
nights and these images were reduced in the same manner as above. Three to six bright, stellar 
objects were used to calibrate the fields observed under nonphotometric conditions. 

3.2. Standards 

The standard stars were processed in an identical manner to that outlined above, though with 
two exceptions. First, the cryocooler pattern was not subtracted from these frames. Tests of the 
images with artificial stars showed that this pattern did not affect the photometry of bright stars. 



*IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 

^Available at http://iraf.noao.edu/iraf/web/contrib.html 
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Secondly, the five to six observations of each standard were photometered individuahy, rather than 
one measurement of a shifted and added frame. The standard stars were measured with PHOT 
using a 15" radius aperture and a mean sky measured in an annulus extending from 25" to 35" 
radius. 

To convert the instrumental magnitudes to the CTIO/CIT system, the fohowing equations 
were solved: 

J = j + k+h (x-i.3)+i2 (J-K-1.4) 

H = h + hQ + hi{X- 1.3) + h2{J -K - 1.4) (1) 
K = k + kQ + kx{X - 1.3) +k2{J-K- 1.4) 

where j, h, k are the instrumental magnitudes, jo, ho, arc the photometric zeropoints, ji, hi, ki 
are the airmass coefficients, and j'l, h'l, k'l arc the color terms. X = 1.3 is the mean airmass of the 
fields observed in this survey and J — K = lA \s approximately the mean J — K color of galaxies 
in the range = 16 — 19 reported by Saracco et al. (1999). While the photometric zeropoints 
and airmass coefficients may vary from night to night, the color terms should not vary so long 
as the filters remain unchanged. To take advantage of the observations of red stars over several 
nights, I therefore performed a multidimensional fit (e.g. Gould 1995) to all of the nights of a given 
observing run simultaneously. For each filter, this allowed the zeropoint and airmass coefficient to 
vary but not the color term. Over a given observing run, the photometric zeropoint varied by a 
few percent. The airmass coefficient was more variable, with large changes generally temporally 
correlated with changes in the weather (as discussed by Progel 1998). If standards at large airmass 
were not observed on a particular night, the average coefficient of the temporally adjacent nights 
on that observing run was used. Most of the fields were observed at X < 1.5 (and all at less than 
two) so uncertainties in the airmass coefficients translate to less than one percent uncertainties in 
the photometry. All of the photometric solution coefficients are listed in Table 2, along with the 
root-mean-square variation in these solutions. 

As mentioned above, there were several changes to the filters in TIFKAM over the course of 
this project. The J filter used in 1997 October had considerable structure in the wings of the 
point spread function (PSF) and therefore it was replaced before the 1998 April observing run. 
The new J filter had a red leak such that the measured sky brightness was Sj !=a 13 mag arcsec"^, 
compared to measurements of Sj ~ 15.5 — 16 mag arcsec"^ on the other observing runs. After 
1998 April this J filter was only used in conjunction with a piece of PK50 glass in the prefilter 
wheel, which serves as a red-blocking filter. PK50 is opaque longward of 2.7/im, but transparent at 
shorter wavelengths. A piece of PK50 was also added to the same filter cell as the H filter between 
the 1998 April and 1998 September observing runs due to the possible presence of a slight red leak 
in this filter. Because of the changes in the J and H filters, the photometric solution has three 
separate color terms for the J filter and two separate color terms for the H filter. The K filter was 
unchanged over the course of this observing program and therefore only one color term was sought 
in this equation. 
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3.3. Object Catalogs 



I cataloged the positions of all objects in these fields with the SExtractor package (Bertin 
&; Arnouts 1996) and the default object detection filter. This choice was motivated by the data 
reduction procedure, which creates images with "dead" regions, corresponding to areas with low 
exposure times and higher sky noise in the final shift-and-add step. The weight map option in 
SExtractor is an efficient way to insure detections on only the useful part of each image file. In 
addition, the noise varies somewhat over our images, and specifically the northern half of the array 
is on average noisier than the southern half. I compared SExtractor to FOCAS (Jarvis & Tyson 
1981; Valdes 1982) and found that SExtractor did a better job of detecting all of the objects in the 
fields without introducing significant spurious detections. 



The limiting magnitude of a survey is typically expressed as the apparent magnitude where 
the probability of detecting an object is equal to some percentage, commonly 50% or 90%. To 
obtain an accurate census of all objects to such limiting magnitudes, one can use an estimate of 
the incompleteness vs. apparent magnitude to correct the observed number of objects to the true 
value. This simple characterization of the detection efficiency is sufficient for the study of objects 
that all have the same size and surface brightness profile and is simplest to apply to point sources 
such as stars or quasars. I discuss the related issue of star-galaxy separation below in §6. 

In studies of galaxies, more compact and higher surface brightness objects are generally easier 
to detect. Not taking the relative "visibility" (Disney 1976; Disney & Phillipps 1983; Phillipps et 
al. 1990; Davies 1990) of different types of objects into account can lead to erroneous conclusions 
about the intrinsic distribution of observed galaxy properties. The detection efficiency of a survey 
should therefore be characterized as not just a function of apparent magnitude, but also as a 
function of surface brightness profile and angular size. That is, surveys should take into account the 
multivariate nature of the galaxy population, rather than collapsing variations in profile shape, size, 
and surface brightness into just a dependence on integrated brightness. Bershady et al. (1998), for 
example, separate objects into stars, small, and large galaxies and calculate the detection efficiency 
of each class separately. 

The surface brightness profile of galaxies is commonly parameterized by either an exponential 
profile 
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Detection Efficiencies 




(2) 



or a de Vaucouleurs (r^/^) profile 




(3) 
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where S(r) is the surface brightness at radius r and r/j is the scale radius for each profile that 
encloses half of the total light. Most galaxies can be fit by either one of these profiles or, for disk 
galaxies with a bulge component, a superposition of both. These paramcterizations are useful for 
examining how the detection efficiency varies as a function of galaxy size and profile shape. 

The sensitivity of each of the fields in this survey is different due to variations in the atmospheric 
conditions and instrument setup. The likelihood of detecting an object of a given brightness, size, 
and profile will vary from field to field and can not be accurately represented by a single number 
for the entire survey, even for stellar objects. To characterize the detection limits of this survey, I 
added artificial objects produced with the IRAF ARTDATA package to each of the fields in each 
filter and measured the fraction of recovered objects with SExtractor. In addition to stellar profiles, 
I added galaxies with both exponential and r^/^ profiles with r/j equal to 0.25", 0.5", 0.75", and 1" 
over a range of apparent magnitude. Studies of the HST Medium Deep Survey (Roche et al. 1997), 
the NICMOS parallel survey (Teplitz et al. 1998; Yan et al. 1998), and HST imaging of the CFRS 
and LDSS galaxies (Lilly et al. 1998) have shown that nearly all galaxies to the apparent magnitude 
limits of this survey are well fit by either exponential or r^^^ profiles and most will have rh < 1" and 
exponential profiles. An alternative to using artificial objects to measure the detection efficiency is 
to extract bright objects from the survey fields, artificially dim them, place them at random back 
into the original images, and attempt to recover them. However, as the median galaxy size declines 
at fainter apparent magnitudes, this may lead to an underestimate of the detection efficiency. 

I convolved the model galaxy profiles used in this analysis with a Moffat (1969) profile as this 
profile is a much better fit to the observed profiles of bright stars than a Gaussian profile. A Moffat 
profile has the functional form 

I{r)=\l + {r/afY . (4) 



After fitting Moffat profiles to bright stars in a number of fields and filters that spanned the range 
of PSF sizes, I fixed /? to be 2.5. The parameter /3 determines the strength of the power-law tail of 
the radial intensity distribution, while a parameterizes the width of the profile peak and is similar 
to C7 in a Gaussian function. This procedure set the convolution kernel for the model galaxy profiles 
in ARTDATA for each image and takes the variation in PSF from image to image into account 
in the determination of the detection efficiency. The functional form of the Moffat profile used by 
the ARTDATA package is slightly different from equation 4, but MOFFAT profiles produced by 
ARTDATA are in reasonable agreement with this standard form. The default dynamic range of the 
convolution kernel in ARTDATA for galaxy profiles (ARTDATA.DYNRANG) is only ten, which 
causes an artificial truncation in galaxy profiles. I changed the value of this parameter to 10000, 
which matches the dynamic range of the PSF computation. While this is more computationally 
intensive, it is not prohibitive. 

In Figure 1 I show the relative detection efficiency for these two profiles at a range of in 
a field of fixed sky brightness. The sky noise and magnitude zeropoint in this figure is from one 
of the K band frames, although the relative detection efficiencies for different sizes, profiles, and 
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Fig. 1. — Variation in the detection elTicicncy vs. apparent magnitude and source size. The upper left panel shows 
the detection efficiency vs. apparent magnitude for a stellar profile (PSF FWHM = 1" and Moffat (3 = 2.5) and four 
exponential disks with intrinsic half-light radii of ru = 0.25", 0.5", 0.75", and 1". The exponential disks have all been 
convolved with the 1" PSF. The lower loft panel is the same except that the intrinsic galaxy surface brightness profile 
is a de Vaucouleur model. The two right panels are from the same image, except with the PSF FWHM = 1.5". The 
sky noise and magnitude zeropoint in this figure is from one of the K band frames, although the relative detection 
efficiencies for different sizes, profiles, and seeing are independent of the actual filter used. 



seeing are independent of tlie actual filter used. Tliroughout this paper I will use Vh to parameterize 
the angular half-light radius, rather than associate it with a physical scale. The top two panels 
of Figure 1 show an exponential disk model and the bottom two panels show a r^/^ profile. The 
left panels correspond to a seeing FWHM of 1" and the right panels correspond to 1.5" seeing. 
This figure shows that the detection efficiency for stellar objects is 0.5 - 1 mag fainter than the 
detection efficiency for marginally resolved galaxies and that the detection efficiency is similar for 
exponential disks and r^/^ profiles with the same r^. As angular size and surface brightness vary 
with redshift, the detection efficiency for different classes of objects can be even more important 
for selecting galaxies for a redshift survey that will be used to compute the luminosity function 
(Dalcanton 1998a; Petrosian 1998). 

Figure 2 illustrates the variation in detection efficiency from field to field of this survey as a 
function of object size. The histogram represents the number of square arcminutes with a given 
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Fig. 2. — Histogram of survey area with a given 90% limiting K magnitude. Each 0.25 magnitude bin shows the 
number of square arcminutes which are 90% complete to that K apparent magnitude for the given object size. The 
top panel shows the distribution for stellar profiles, while the remaining panels show the distribution for exponential 
disks with half-light radii of = 0.25", 0.5", 0.75", and 1". 

90% completeness limit vs. K magnitude for the 185 square arcminutes of this survey. The top 
panel shows the distribution for stellar profiles, while the lower panels show exponential disks with 
Th = 0.25", 0.5", 0.75", and 1". The average, area-weighted 50% completeness limits of the survey 
are J = 20.5, H = 19.5, and if = 18 mag for exponential disks with = 0.75". For stellar objects 
these limits are J = 21, H = 20, and K = 18.5 mag. 

5. Photometry 

The unknown intrinsic morphology of galaxies makes galaxy photometry considerably more 
challenging than measuring the apparent magnitude of a star. In this paper, the goal is to measure 
the total, integrated brightness of galaxies to compare to models of galaxy surface density as a 
function of integrated brightness. Most of the galaxies in this sample have small angular sizes, 
where small in this context corresponds to rh less than or equal to the seeing FWHM. Due to the 
finite signal-to-noise ratio and the lack of numerous bright stars in each field to exactly determine 
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the PSF, the true morphology of the galaxies can not be accurately deconvolved. Photometric 
techniques for marginally-resolved galaxies include aperture magnitudes, isophotal magnitudes, 
metric magnitudes, and various slight modifications of these techniques. 

Aperture magnitudes are the simplest to define as they just involve computing the total flux 
within some given size aperture. The problem with aperture magnitudes, however, is that they will 
always exclude some fraction of a galaxy's light unless the aperture is made very large. While using 
a very large aperture is a way of avoiding systematic errors in galaxy photometry, it is impractical to 
apply to observations within several magnitudes of the survey limit due to nonnegligible sky noise. 
A pragmatic approach is therefore to measure the brightness of the galaxy within some small 
aperture and then to correct this aperture magnitude to compensate for the lost light, usually 
based on some measure of the galaxy size (Glazebrook et al. 1994; Soifer et al. 1994; Cowie et al. 
1994; Bershady et al. 1998; Saracco et al. 1999). For any size aperture, the fraction of lost light 
will depend on both the scale size and surface brightness profile of the galaxy. These quantities are 
difficult to measure with great accuracy for galaxies near their detection limit. 

An isophotal magnitude is measured by summing all of a galaxy's flux out to some limit- 
ing isophote, usually the limiting isophotal magnitude of the survey. A variation of the isophotal 
magnitude called the "TOTAL" magnitude is more commonly used (McLeod et al. 1995; Hall et 
al. 1998; Minezaki et al. 1998b) and is one of the quantities calculated by FOG AS. The FOG AS 
TOTAL magnitude is defined as the sum of the flux within twice the isophotal area, rather than 
just the isophotal area. The isophotal magnitude is well-known to be biased towards underesti- 
mating the flux of all but very bright and high surface-brightness galaxies. Using the FOCAS 
TOTAL magnitude decreases the tendency of isophotal magnitudes to underestimate galaxy size 
and integrated brightness by effectively increasing the aperture size, but it is fundamentally still 
an isophotal magnitude and therefore underestimates the integrated light of galaxies (Dalcanton 
1998b). 

Metric magnitudes appear to offer a physically motivated solution to the shortcomings of these 
two techniques. The most common metric magnitude is measured within the radius introduced by 
Petrosian (1976). The Petrosian magnitude is defined to be the flux interior to radius Vp, where 
the surface brightness at Vp is equal to some fraction of the average surface brightness interior to 
rp. The great attraction of using a Petrosian magnitude is that it is independent of the spacetime 
geometry by virtue of being the ratio of two surface brightnesses. Vp therefore corresponds to the 
same physical scale in two galaxies with the same surface-brightness profile at different redshifts. 
The corresponding integrated brightness within that radius will correspond to the same fraction 
of the total light (any systematic underestimate will be identical for the two galaxies). However, 
the fraction of the light missed by a Petrosian magnitude for two galaxies at different redshifts 
will only be the same if two conditions hold: the surface brightness profiles of the galaxies are the 
same and neither of the galaxies are small enough that they are significantly affected by the PSF. 
Thus while the Petrosian magnitude is ideal for cosmological tests (Sandagc &: Perelmuter 1990; 
Petrosian 1998), it may not offer the ideal choice for measuring the integrated light of galaxies with 
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a diverse range of surface brightness profiles and sizes comparable to the PSF of the observations. 

The main goal of the magnitude measurements in this paper is to measure the number of 
galaxies per unit magnitude. These measurements, in conjunction with theoretical predictions, can 
then be used to constrain models of galaxy formation and evolution. The output of the theoretical 
models to which observed number counts are compared is the sum of galaxies in a given apparent 
magnitude bin. To compare observational data to these theoretical predictions therefore requires 
an unbiased measurement of the number of galaxies at a given integrated brightness, which in turn 
requires a good means of estimating the integrated brightness of individual galaxies. 

I evaluated different techniques for measuring the integrated brightness of the galaxies in 
the sample including aperture magnitudes measured with PHOT in IRAF, TOTAL magnitudes 
measured with FOCAS, and "BEST" magnitudes measured by SExtractor. The SExtractor BEST 
magnitude is defined to be either a Kron magnitude (Kron 1980), which is measured in an elliptical 
aperture whose size is determined by the object's profile, or an isophotal magnitTidc in very crowded 
regions. To evaluate the accuracy of these photometry techniques I used ARTDATA in IRAF (as 
described in §4) to add artificial galaxies to the images. The metric to evaluate the accuracy of 
each of these techniques is the size of the difference between the input magnitude and the measured 
magnitude as a function of rh for exponential and r^/^ profiles. I performed these tests for different 
images to sample the range of seeing in the data and at a range in apparent magnitude from the 50% 
completeness limit for stellar profiles to four magnitudes brighter. This magnitude range includes 
more than 80% of the sources in the sample. The investigation showed that IRAF PHOT aperture 
magnitudes, FOCAS TOTAL magnitudes, and SExtractor BEST magnitudes all underestimate 
the integrated brightness of galaxies. The FOCAS and SExtractor magnitudes were taken directly 
from the output of these two packages, while I added a stellar-profile aperture correction, based on 
measurements of several stellar objects in the field, to the aperture magnitudes measured by IRAF 
PHOT. 

Figure 3 shows the result of one of these experiments for an image with 1.25" FWHM. The 
measured magnitude is either the IRAF PHOT aperture magnitude (left panels), FOCAS TOTAL 
magnitude {middle panels), or SExtractor BEST magnitude {right panels) and the true magnitude 
is the input magnitude assigned to the artificial object. As in Figure 1, the sky noise and apparent 
magnitude for this simulation are based on one of the K band frames; however, the reliablility 
of these different photometric techniques relative to the completeness limit are independent of 
the filter. This figure shows the results for stellar profiles {top panels), = 0.5" exponential 
disks {middle panels), and = I" {bottom panels). The data points are the mean offset between 
the measured and true magnitudes from ~ 200 artificial objects per apparent magnitude and the 
errorbars are the la dispersion in these measurements at that magnitude. At a fixed input galaxy 
brightness more light is missed in galaxies with larger rh, while at fixed rh, more light is missed 
in an r^^'^ profile than in an exponential profile. For fainter apparent magnitudes, the size of the 
brightness underestimate increases for the FOCAS TOTAL and SExtractor BEST magnitudes at 
fixed rh, but remained constant for the aperture magnitude. The reason for the systematic increase 
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in the underestimate for the FOCAS TOTAL magnitudes is that the isophotal area of a galaxy of 
fixed size decreases as the integrated brightness decreases. Though the FOCAS TOTAL magnitude 
doubles the isophotal area, this still leads to an underestimate of the true size of the galaxy for 
the photometric measurement. The SExtractor BEST magnitude suffers from the same systematic 
increase in error as the size of the elliptical aperture is set by a fit to the observed galaxy profile. 
This offset is also due to the fact that SExtractor BEST magnitudes assume the PSF is Gaussian 
when computing the elliptical aperture for photometry and a Gaussian profile is more centrally 
concentrated than a MOFFAT profile (Dalcanton 1998b). At a higher galaxy surface density than 
probed in this survey, the SExtractor BEST magnitudes will be isophotal magnitudes, rather than 
Kron magnitudes, and the brightness underestimate may be more severe. For fainter galaxies at 
fixed size and profile shape, a smaller fraction will be above the noise and the tendency is for a 
detected object to appear more compact, which will further shrink the aperture size used in the 
photometry. While I have not evaluated Petrosian magnitudes here, Dalcanton (1998a) included 
them in her study of systematic biases in measuring the integrated light of galaxies for computing 
the luminosity function. She found that Petrosian magnitudes underestimate the true brightness 
of galaxies to a lesser extent than aperture or isophotal magnitudes (see also Dalcanton 1998b), 
although they too are subject to large errors for objects near the detection limit. This is because the 
surface brightness profile becomes increasingly noisy. In addition, as stated above, the Petrosian 
radius is difficult to determine for marginally resolved objects. 

Saracco et al. (1999) also compare aperture magnitudes to isophotal magnitudes and SExtrac- 
tor AUTO magnitudes, which are similar to SExtractor BEST magnitudes. Their figure 1 shows 
a difference of several tenths of a magnitude between photometric measurements with these three 
techniques of the galaxies in their fields. These differences increase to on order one magnitude 
for the last two magnitudes below their S/N = 5 detection limit for compact objects. They find 
that both isophotal and SExtractor AUTO magnitudes underestimate the flux of faint objects and 
therefore chose to use aperture magnitudes. Vaisanen et al. (2000) have also evaluated the relative 
merits of aperture, isophotal, and SExtractor BEST magnitudes. They concluded that SExtractor 
BEST magnitudes do not underestimate the integrated brightness of galaxies near the detection 
limit, in contrast with the simulations presented here and the results of Dalcanton (1998b). 

My analysis has shown that only aperture magnitudes avoid errors in the photometry that 
are a function of integrated brightness, and I have chosen to use them. The brightness underes- 
timate of aperture photometry will still depend on the intrinsic size and profile shape, but avoids 
the additional systematic uncertainty due to the dependence on integrated brightness inherent in 
SExtractor BEST and FOCAS TOTAL magnitudes. Another consideration in the decision to use 
aperture magnitudes is the unique nature of these data. I have 21 fields observed in three filters, 
each with a different PSF shape. Aperture photometry is the best way to account for the hetero- 
geneous nature of the dataset. The size of the error in the SExtractor BEST and FOCAS TOTAL 
magnitudes will vary from field to field with the seeing. While aperture magnitudes will suffer from 
the same problem, they can be relatively easily scaled to take seeing variations into account. In 
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Fig. 3. — Reliability of the photometric methods as a function of image size and integrated brightness. The left 
panels show the difference between the measured and true magnitudes for aperture photometry with IRAF PHOT. 
The top panel shows the results for stellar objects, the middle and bottom panels for exponential disks with rh = 0.5" 
and Vh = 1", respectively. The middle panels show the same experiment with FOCAS TOTAL magnitudes and the 
right panels with SExtractor BEST magnitudes. The dashed, vertical line is the 90% detection limit for the objects 
in each panel. As in Figure 1, the sky noise and magnitude zeropoint were taken from one of the K band frames; 
however, these results are only dependent on the apparent magnitude relative to the magnitude of the 90% detection 
limit. 



addition, the aperture magnitudes can be easily measured through a sufficiently large aperture to 
minimize the size of the galaxy aperture correction and therefore the associated uncertainty. 

Given the range of galaxy sizes and typical seeing, aperture magnitudes with a radius 1.5 
times the FWHM of the seeing are a good compromise between minimizing the size of the galaxy 
aperture correction and maximizing the galaxy signal in the aperture. At z = 0.5, 1" corresponds 
to a physical size of 4.3 kpc (for Qm = 0.3, Qa = 0.7) and therefore these typical aperture sizes 
of 1.5" - 2.5" include most of the galaxies' light. The galaxy aperture correction is the difference 
between the true integrated brightness of the galaxy and the measured aperture magnitude plus the 
stellar aperture correction. The stellar aperture correction only depends on the shape of the PSF 
and, because the aperture is scaled to the PSF size, this correction was always measured to be ~ 0.1 
mag. The galaxy aperture correction depends on the typical galaxy size and was measured with 
the same code developed to measure the detection efficiency for objects of different size, surface 
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brightness profile, and integrated brightness. The distribution of size and surface brightness profile 
as a function of integrated brightness have not been completely determined in the NIR. However, 
a range of studies of galaxy morphology at visible (Small et al. 1995; Roche et al. 1997, 1998) 
and NIR (Yan et al. 1998; Bershady et al. 1998; Saracco et al. 1999) wavelengths, combined with 
color information (Thompson et al. 1999) shows that most of the faint galaxies in this survey will 
have exponential profiles with typical half-light radii of 0.5" - 0.75". At the magnitude limits of 
this sTirvey the average galaxy size is ^ 0.6" in these three filters. Any bulge component for 
these galaxies, though more prominent at NIR wavelengths, will fall inside the aperture. For this 
range of sizes of galaxies, I computed the galaxy aperture correction for each frame as a function 
of integrated brightness. The size of this correction is always in the range 0.1 - 0.2 mag. The 
correction is generally smaller for fields with poorer seeing as the aperture size is larger in these 
fields, although the detection limit for a given galaxy size and profile shape is lower for these fields 
as well. 



6. Star/Galaxy Separation 

The calculation of number counts of galaxies at faint magnitudes must deal with stellar con- 
tamination, which affects faint galaxy surveys in two ways. First, the typical galaxy size decreases 
at fainter apparent magnitudes and therefore galaxies are progressively more likely to be morpho- 
logically indistinguishable from stars. Second, stars (and compact galaxies) are the easiest objects 
to detect and accurately photometer near the survey limits. This is illustrated in Figure 1, which 
shows that the limiting magnitude at fixed detection efficiency is on order one magnitude brighter 
for exponential disks or r^/^ profiles with rh = 1" than for stars. Multicolor separation is one way 
to remove stellar contamination from faint galaxy samples. This process is most effective with small 
photometric errors, although it can be effective with larger errors with a longer wavelength baseline 
(e.g. visible-NIR colors, Gardner 1995). The removal of stellar contamination from this survey is 
particularly important as stars contribute on order 20 - 30% to even the faintest magnitude bins 
of this survey. In even deeper surveys stellar contamination is less of a problem as galaxies begin 
to significantly outnumber stars. 

The stellar content of the DMS fields has already been cataloged as part of the search for 
quasars (Hall et al. 1996a; Osmer et al. 1998). The selection criterion for identifying stellar objects 
was that they must be indistinguishable from the PSF in at least three of the six CCD filters (Hall 
et al. 1996a). I used the stellar catalog from Osmer et al. (1998) to identify and filter out stars 
from the NIR catalogs. Given that the CCD stellar catalog extends several magnitudes fainter than 
this survey, all but the very reddest stars should be included in the stellar catalog. In contrast, 
the CCD stellar catalog includes some contamination from compact galaxies near its magnitude 
limit. From the stellar model counts presented in Minezaki et al. (1998a) and galaxy counts in the 
literature, I expect stars comprise approximately 10 - 20 % of all objects at the limiting magnitude 
of this survey. Any missed fraction will therefore be a minor contamination to the galaxy number 
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counts. 

7. Number Counts 

The calculation of NIR galaxy number counts in the DMS is complicated by the variation in 

detection efficiency from field to field. To combine all of the fields I have used a technique analogous 
to the one outlined in the appendix of Bershady ct al. (1998). However, the same variations in 
seeing that contribute to the different field sensitivities make it difficult to reliably bin these objects 
into different size classes. I have therefore not separated the counts into different sizes (and profile 
types), effectively collapsing the galaxy distribution from a multivariate distribution (size, profile 
shape, and integrated brightness) to a fTinction of integrated brightness alone. As most number- 
magnitude models for galaxies only predict surface density vs. integrated brightness (though see 
Im et al. 1995), this is the most useful format to compare the data with theoretical predictions. 

For each apparent magnitude bin I used the typical galaxy size and galaxy aperture correction 
discussed in §5 to correct the measured brightness of each galaxy. All of the galaxies in each field 
were then placed into 0.5 mag bins in apparent magnitude. The detection efficiency for each bin is 
the average of the detection efficiency vs. magnitude measured in 0.1 mag increments and weighted 
by the expected slope of the number-magnitude relation a ~ 0.5, where logiV (x am. In practice, 
the detection efficiency drops off significantly faster than the number-magnitude relation increases 
and so the detection efficiency of a magnitude bin is relatively insensitive to a. If the detection 
efficiency of the faint end of a bin was less than 50%, this bin was not included in the calculation 
of the number-magnitude relation. The detection efficiencies described in §4 multiplied by the 
field size yields the effective area of each field. I then summed the total differential counts and 
effective area of each field per 0.5 magnitude bin to compute the differential number counts per 
unit magnitude per square degree. The number counts for J, H, and K are shown in Figures 4, 
5, and 6, respectively. The raw and corrected counts, along with the effective area, are listed in 
Tables 3, 4, and 5. The raw counts are the total counts in 0.5 mag bins with only the stellar (~ 0.1 
mag) aperture correction applied to the photometry. The corrected counts are the counts per 0.5 
mag bin after application of the galaxy aperture correction. The errorbars on the number counts 
are la upper and lower confidence intervals calculated using the formulae in Gehrels (1986). 

The J number counts shown in Figure 4 agree well with the number counts of Teplitz et al. 
(1999), Vaisanen et al. (2000) and the (small+large) counts of Bershady et al. (1998), though they 
are somewhat higher than the Saracco ct al. (1999) counts. The Vaisanen ct al. (2000) survey 
overlaps the most with this survey, although they use SExtractor BEST magnitudes rather than 
aperture magnitudes. Their survey covers a substantially larger area than this one, though its faint 
limit is nearly one magnitude brighter. Bershady et al. (1998) and Saracco et al. (1999) extend 
much fainter than this survey and their brightest counts only overlap with the faintest magnitudes 
of this survey and in this regime they have approximately five to ten times fewer galaxies due to 
their smaller area. The Saracco et al. (1999) surface density is probably lower in this range due to 
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Fig. 4. — DifFcrcntial galaxy number counts in the noar-infrarcd J band. The raw number counts have been 
converted into units of mag~^ deg~^ (filled circles) and the errorbars are la confidence limits. These data are listed 
in Table 3. Also shown are number counts from the literature and four models from Gardner (1998). 

small number statistics as they also use the Persson et al. (1998) standards, aperture magnitudes, 
and the agreement appears to improve at fainter magnitudes where they have a larger sample. All 
of the J counts in the range 16 - 20.5 mag have essentially the same slope as these data, a ~ 0.54. 

The slope of the H band counts in Figure 5 over the 16.5 - 19 mag range is a = 0.47. These 
counts overlap with the brightest magnitude bin from Yan et al. (1998), but the surface density 
vs. magnitude relation shows an apparent "break" at = 19 - 20 mag. The Yan et al. (1998) 
measurements are based on NICMOS parallel data with Camera 3 over a total area of 8.7 arcmin^ 
and therefore they have relatively few galaxies in the region of overlap. However, the additional 
NICMOS parallel data studied by Teplitz et al. (1998), which includes images from Camera 2 and 
extends several magnitudes brighter, also shows an apparent break at i7 = 19 - 20 mag. Teplitz et 
al. (1998) transformed several deep K galaxy counts studies to the NICMOS F160W filter and these 
transformed K counts appear to match the break in the H galaxy number-magnitude relation. This 
break may therefore be the change in slope seen in K band counts at ii' ~ 20 mag. 

Most published NIR number counts to date have been obtained in the K band (Gardner et al. 
1993; Cowie et al. 1994; Glazebrook et al. 1994; Soifer et al. 1994; Djorgovski et al. 1995; McLeod 
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Fig. 5. — Differential galaxy number counts as in Figure 4 but for the near-infrared H band. 

ct al. 1995; Moustakas et al. 1997; Minczaki et al. 1998b; Szokoly et al. 1998; Bcrshady et al. 
1998; Saracco et al. 1999; Vaisanen et al. 2000; McCracken et al. 2000). Gardner et al. (1996) and 
Vaisanen et al. (2000) have the largest magnitude range of overlap with the DMS counts shown 
in Figure 6 and the number-magitude relations from these studies are in good agreement. The 
slope of the K counts from 14 - 18 mag is a = 0.54. At fainter magnitudes there is a great deal 
of dispersion in the surface density measurements, though given the small areas of these surveys 
they are generally consistent with one another. As noted by several authors (e.g. Bershady et al. 
1998; Saracco et al. 1999), this dispersion may also be due to different measurement techniques, 
variations in magnitude systems, and filters, in addition to small number statistics and cosmic 
variance. The most significant differences are likely due to the method of magnitude measurement, 
completeness corrections, and at least at bright magnitudes, star-galaxy separation. Of the faint 
counts presented in Figure 6, Saracco et al. (1999) and Moustakas et al. (1997) apply a uniform 
galaxy aperture correction for all objects, Bershady et al. (1998) compute completeness and a 
galaxy aperture correction as a function of object size, Minezaki et al. (1998b) use FOCAS TOTAL 
magnitudes, and Djorgovski et al. (1995) only use a stellar aperture correction. 

I computed the four models shown in Figures 4-6 with ncmod, a publically available code 
for computing galaxy number counts and colors by Gardner (1998). The models correspond to no 
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evolution (solid lines), passive evolution (dotted lines), merging plus passive evolution (short-dashed 
lines), and dust plus passive evolution (long-dashed lines). For all of the models I have assumed 
CIm = 0.3 and JIa = 0.7, in agreement with results on Type la supernovae (Riess et al. 1998; 
Perlmutter et al. 1999), the Lyman-a forest, COBE-DMR (Phillips et al. 2000), as well as number 
counts in the Hubble Deep Field (Totani & Yoshii 2000). I chose to fix the cosmological parameters 
as the uncertainties in ^Im and J7a are less than the uncertainties in models of galaxy evolution. 

All of these models are based on the K luminosity function and mix of galaxy types from 
Gardner et al. (1997) and galaxy spectral energy distributions from GISSEL96 (Bruzual & Chariot 
1993). To compute these models I also used the filter transmission curves for the TIFKAM filters 
plus a measurement of the atmospheric transmission. The no evolution model is a pure fc-correction 
model. In the passive evolution model, all galaxies form at z = 15 except for the latest galaxy type, 
which have constant star formation and are always 1 Gyr old. The merging model has the same 
galaxy evolution as the passive evolution mode, but it also includes number evolution in the form 
0* oc (1 + z)^ with 7] = 1.5 and conservation of luminosity density. This model is based on the 
parameterization of Rocca-Volmerange &: Guiderdoni (1990) and it produces similar results to the 
Broadhurst et al. (1992) merging model. This value of r/ is higher than the constraint rj < 1 derived 
by Totani Sz Yoshii (2000) in their A model for the HDF, but is illustrative of the effect of adding 
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merging to passive evolution. Finally, the dust model also contains the same galaxy evolution 
parameterization as the passive evolution model, but with the addition of dust in the plane of the 
galaxies. This model, based on the work of Bruzual ct al. (1988) and Wang (1991) was introduced 
to number count models to correct for their tendency to overproduce the UV flux in galaxies, 
though it has a minimal effect on the NIR number-magnitude relation. 

The J, H, and K observations presented here all agree well with the passive evolution model. 
They are less consistent with the number density evolution we assumed in the merging model, and 
clearly have a higher surface density of objects than the no evolution prediction. The observed J 
counts are all higher than the models, but this may be due to the extrapolation of the K galaxy 
luminosity function to the J band. As noted above, the HST counts for H > 20 mag (Teplitz et al. 
1998; Yan ct al. 1998) have systematically lower surface density. This difference does not appear 
to be due to differences in the transmission profiles of the two filters. When I calculated the model 
galaxy counts with the F160W transmission profile, rather than the H filter, the F160W galaxy 
count models were offset by less than 0.1 mag fainter than the ground-based H model counts, 
whereas an offset of 0.3 - 0.4 mag is needed to explain the difference in surface density. 

8. Mean Colors 

Color as a function of apparent magnitude can provide an additional constraint on galaxy 
evolution models, particularly when severals colors are available. I have used the Gardner (1998) 
model to compute the median J ~ H color vs. J magnitude and the median J — K color vs. K 
magnitude for the four models discussed above. Figure 7 shows these two color-magnitude relations 
for all of the galaxies detected in these filters. For the magnitude range of this survey, all of these 
models predict similar median colors and they agree well with the data and J — K colors reported 
by Bershady et al. (1998) and Saracco et al. (1999) in the region of overlap. 

To these apparent magnitude limits, NIR colors do not provide a strong constraint on galaxy 
evolution as most galaxies are low redshift {z < 1) objects. In this redshift range, the JHK filters 

all sample the old stellar population, which has a relatively flat spectral energy distribution inde- 
pendent of galaxy type. The overlap of the four models and the lack of change in the median color 
with apparent magnitude reflects the similarity of these galaxies at NIR wavelengths. Bershady et 
al. (1998) and Saracco et al. (1999) also report J — K colors vs. K and these results agree well with 
their colors. 

9. Summary 

I have presented NIR JHK galaxy counts and colors and described in detail the detection 
efficiencies and selection effects that affect this survey. I have studied several popular methods 
for measuring the integrated brightness of galaxies as a function of apparent magnitude and found 
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Fig. 7. — Mean near- infrared colors vs. magnitude for the galaxies in this sample. The data are represented by 
points with errorbars. Only every tenth point has an associated errorbar for clarity. The curves are the mean colors 
vs. magnitude for the four models shown in Figures 4 6 and discussed in §7. The straight, dashed lines arc the 
average 50% limiting magnitudes and colors for an exponential disk with rn = 0.75". The scatter of the data points 
about the model lines is primarily due to variance in the typical galaxy colors, rather than photometric uncertainty. 

that only aperture magnitudes miss a consistent fraction of a galaxy's light as a function of ap- 
parent magnitude for a range of galaxy sizes and profile shapes. This survey is the largest area 
observed to these depths in all three NIR filters that has been published to date. The J, H, and 
K number-magnitude relations and colors are in good agreement with a simple model of passive 
galaxy evolution with at most a small amount (r/ < 1.5) of merging in an Q^m = 0.3, i^A = 0.7 
cosmology. 

Current model predictions for the number-magnitude relation at NIR wavelengths are very 
similar to one another and larger surveys at the faintest magnitudes are needed to discriminate 
between different models with this diagnostic tool. The models shown in Figures 4-6 suggest that 
number counts to H ^ 2A mag or i^' 23 mag are necessary to distinguish between the passive 
evolution and merging model. High angular resolution imaging to measure the mix of galaxy sizes 
and morphological types as a function of brightness, as well as color information, could prove to be 
a more efficient means to discriminate between galaxy formation and evolution models with deep, 
NIR imaging surveys. 
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Table 1. Log of Observations 
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0. 1 L 


CiQCiA 1 1 


non/i 1 
yoU4iO 


980410 


1.5 


i.i 


1 




y.i / 


yoU4i4 


noH/i 1 /I c 
yoU4i4 


08041 4^ 


1 3 

J- - fJ 


1 
I.Z 


i.D 




y.zo 


yoiooo 


noi nno 
yolOUo 


981003 


1.6 


1 ^? 

1.0 


1 
I.Z 




o.oi 


OQi nno 


noi nno 
yoiUUz 


981002 


1.3 


1 
I.z 


1 

1.0 




y.y4 


no 1 nno 
yolOOo 


noi nno 

yoiooz 


981002 


1.3 


1 

1.0 


1 cc 
1.5 


21EC 


9.31 


981001 


981001 


981002 


1.4 


1.2 


1.5 


01WC150W 


9.86 


981003 


981003 


981003 


1.4 


1.5 


1.5 


21WC150W 


9.30 


981004 


981004 


981004 


1.7 


1.5 


1.4 


14NC150W 


9.04 


990516 


990516 


990516 


1.0 


1.0 


1.0 


14NC150E 


9.36 


990517 


990517 


990517 


1.6 


1.0 


1.2 


14NC300W 


9.90 


990519 


990519 


990519 


1.2 


1.1 


1.2 


14NC300E 


9.54 


990520 


990520 


990520 


0.9 


0.9 


1.0 


17SC 


9.70 


990515 


990926 


990927 


1.1 


1.1 


1.1 


21WC150E 


9.80 


990928 


990928 


990928 


1.3 


1.3 


1.3 


22EC150W 


10.19 


990929 


990929 


990929 


1.7 


1.8 


1.9 



^Area that does not overlap 01WC150W 

''Calibration obtained 981003 
'^Calibration obtained 990515 

Note. — Log of the observations of the 21 subfields of the DMS. Column 1 lists the 
field identifier, which indicates either the center of one of the fields listed in Hall et al. 
(1996a) or the offset from the center of one of these fields. For example, 21WC refers to 
the center of Field 21W and 01EC60S refers to a field 60" south of Field OlE. The two 
exceptions to this naming convention are CFl and CF3, which are defined in Hall et al. 
(1998). Column 2 has the area in square arcminutes of each subfield, while columns 3 - 
5 the UT date of the observations in J, H, and K, respectively. Columns 6-8 list the 
FWHM in arcseconds of the PSF in the final, combined frame for J, H, and K. 
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Table 2. Photometric Solutions 



Date jo jl j2 o-(J) ho hi /12 cr(i?) fco h h a{K) 



971011 


22, 


,204 


-0. 


,040 


0, 


,160 


0, 


.046 


21. 


,961 


0.000 


0, 


.060 


0, 


,074 


21 


.245 


-0. 


,060 


0, 


,000 


0, 


.127 


971014 


22. 


,303 


-0. 


,047 


0, 


.160 


0, 


.036 


22. 


,040 


-0.006 


0, 


.060 


0, 


.048 


21, 


.290 


-0. 


,060 


0, 


.000 


0, 


.048 


971015 


22. 


.351 


-0. 


,038 


0, 


.160 


0, 


.048 


22. 


,076 


-0.021 


0, 


.060 


0, 


.029 


21, 


.354 


-0. 


,053 


0, 


.000 


0, 


.077 


971016 


22. 


.346 


-0. 


,040 


0, 


.160 


0, 


.036 


22. 


,063 


-0.040 


0, 


.060 


0, 


.034 


21, 


.324 


-0. 


,060 


0, 


.000 


0, 


.064 


971017 


22. 


.330 


-0. 


,040 


0, 


.160 


0, 


.031 


22. 


,070 


0.000 


0, 


.060 


0, 


.035 


21, 


.320 


-0. 


,060 


0, 


.000 


0, 


.039 


980410 


22, 


,425 


-0. 


,113 


0, 


.060 


0, 


.034 


22. 


,033 


-0.058 


0, 


.060 


0, 


,036 


21 


.316 


-0. 


,050 


0, 


.000 


0, 


.035 


980411 


22, 


,409 


-0. 


,206 


0, 


.060 


0, 


.022 


22, 


,047 


-0.158 


0, 


.060 


0, 


.028 


21 


.311 


-0. 


,134 


0, 


.000 


0, 


.048 


981001 


22, 


,293 


-0. 


,060 


0, 


.060 


0, 


.032 


22. 


,075 


-0.037 


0, 


.000 


0, 


.021 


21 


.372 


-0. 


,065 


0, 


.000 


0, 


.032 


981002 


22. 


.336 


-0. 


,070 


0, 


.060 


0, 


.019 


22. 


,147 


0.000 


0, 


.000 


0, 


.035 


21, 


.485 


-0. 


,052 


0, 


.000 


0, 


.017 


981003 


22. 


.355 


-0. 


,081 


0, 


.060 


0, 


.032 


22. 


,164 


-0.026 


0, 


.000 


0, 


.020 


21, 


.485 


-0. 


,091 


0, 


.000 


0, 


.033 


981004 


22. 


.332 


-0. 


,102 


0, 


.060 


0, 


.021 


22. 


,165 


-0.058 


0, 


.000 


0, 


.022 


21, 


.470 


-0. 


,113 


0, 


.000 


0, 


.024 


990515 


22, 


,377 


-0. 


,090 


0, 


,060 





.040 


22. 


,182 


0.000 


0, 


.000 


0, 


.017 


21 


.494 


-0. 


,050 


0, 


,000 


0, 


.029 


990516 


22, 


.412 


-0. 


,090 


0, 


,060 


0, 


.032 


22. 


,181 


-0.008 


0, 


.000 


0, 


.025 


21 


.505 


-0. 


,046 


0, 


,000 


0, 


.019 


990517 


22, 


,389 


-0. 


,090 


0, 


.060 


0, 


.023 


22. 


,175 


-0.032 


0, 


.000 


0, 


.034 


21 


.495 


-0. 


,043 


0, 


.000 


0, 


.026 


990519 


22. 


.430 


-0. 


,090 


0, 


.060 


0, 


.023 


22. 


,173 


-0.031 


0, 


.000 


0, 


.029 


21, 


.494 


-0. 


,083 


0, 


.000 


0, 


.036 


990520 


22. 


.425 


-0. 


,088 


0, 


.060 


0, 


.041 


22. 


,149 


-0.007 


0, 


.000 


0, 


.023 


21, 


.472 


-0. 


,091 


0, 


.000 


0, 


.020 


990926 


22, 


,204 


-0. 


,047 


0, 


.060 


0, 


.037 


22. 


,030 


-0.006 


0, 


.000 


0, 


,040 


21 


.342 


-0. 


,100 


0, 


,000 


0, 


.071 


990927 


22, 


,208 


-0. 


,147 


0, 


,060 


0, 


.016 


22. 


,045 


-0.050 


0, 


.000 


0, 


.014 


21 


.363 


-0. 


,108 


0, 


.000 


0, 


.031 


990928 


22, 


,211 


-0. 


,130 


0, 


.060 


0, 


.030 


22. 


,054 


-0.064 


0, 


.000 


0, 


.019 


21 


.392 


-0. 


,105 


0, 


.000 


0, 


.064 


990929 


22, 


,240 


-0. 


,125 


0, 


.060 


0, 


.028 


22. 


,072 


-0.073 


0, 


.000 


0, 


.021 


21 


.381 


-0. 


,110 


0, 


.000 


0, 


.026 



Note. — Photometric solutions for all of the clear nights. Column 1 lists the UT date of the observations, 
while columns 2 - 4, 6 - 8, and 10 - 12 list the coefficients of the photometric solutions defined in §3.2. Columns 
5, 9, and 13 list the rms scatter in the solution for each filter and each night. 
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Table 3. Differential J Number Counts 



Mag 


N 


Area 


N 

^ ' corr 


Nlow 


N 


15.50 


2 


0.05150 


117 


53 


230 


16.00 


2 


0.05150 


78 


28 


180 


16.50 


4 


0.05150 


155 


81 


278 


17.00 


3 


0.05150 


272 


172 


418 


17.50 


28 


0.05150 


1049 


848 


1292 


18.00 


29 


0.05037 


1390 


1156 


1666 


18.50 


64 


0.04453 


3144 


2769 


3566 


19.00 


114 


0.04022 


5868 


5329 


6460 


19.50 


145 


0.03041 


10261 


9440 


11150 


20.00 


138 


0.02356 


12394 


11370 


13507 


20.50 


59 


0.00753 


15137 


13139 


17419 



Note. — Differential J number counts of galaxies. Column 1 
lists the center of each 0.5 magnitude bin, while column 2 lists the 
raw number of galaxies detected in that magnitude range. Column 
3 lists the total effective area of that magnitude range in square 
degrees and column 4 contains the number of galaxies per magnitude 
per square degree after accounting for the aperture correction and 
detection probability. Columns 5 and 6 show the la lower and upper 
confidence limits, respectively, computed as in Gehrels (1986). 
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Tablc 4. Differential H Number Counts 



Mag Nraw Area iVcorr ^low ^upp 



15.00 2 0.05310 75 27 175 

15.50 4 0.05310 151 79 270 

16.00 3 0.05310 113 52 223 

16.50 12 0.05310 791 620 1004 

17.00 28 0.05226 1301 1079 1565 

17.50 59 0.05145 2060 1778 2384 

18.00 97 0.05132 4131 3730 4572 

18.50 130 0.04363 6463 5919 7054 

19.00 122 0.02187 11524 10499 12644 



Note. — Same as Table 3 for the differential H number counts. 



Table 5. Differential K Number Counts 



Mag Nraw Area iVcorr Niow 



14.00 2 0.05310 75 27 175 

14.50 3 0.05310 113 52 223 

15.00 4 0.05310 113 52 223 

15.50 5 0.05310 339 228 494 

16.00 21 0.05248 1029 832 1268 

16.50 50 0.05163 2169 1880 2499 

17.00 59 0.04617 3379 2997 3807 

17.50 81 0.03438 3607 3150 4126 

18.00 50 0.01424 7163 6163 8312 



Note. — 



Same as Table 3 for the differential K number counts. 



